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ABSTRACT 


The  f igure-of-merlt  form  of  the  sonar  equation  causes  a 
bias  in  calculated  broadband-detection  ranges.  The  bias  is 
caused  by  the  necessity  of  treating  the  figure  of  merit  as  if 
it  were  concentrated  at  a single  frequency,  since  only  single- 
frequency— transmission-loss  graphs  are  practically  available. 

A signal-excess  form  of  the  sonar  equation  is  used  in  this  report 
to  remove  the  bias  from  passive,  broadband,  aural,  detection- 
range  calculations. 


INTRODL'CTION 

The  purpose  of  this  paper  is  to  develop  a method  for  calculating 
passive,  broadband,  aural  detection  ranges  that  is  more  reliable  than 
the  present  method.  Broadband  detections  are  important  since  many,  if 
not  most,  detections  are  made  on  broadband  energy. 

In  part  1 of  this  report,  the  method  will  be  developed  and  illus- 
trated with  an  example.  In  Part  2,  more  detailed  examples  will  be  pre- 
sented. 

In  almost  all  present  practice,  a single-frequency  method  is  used 
for  calculating  detection  ranges.  Generally,  all  terms  of  the  sonar 
equation  are  functions  of  frequency.  Adequate  treatment  of  the  broadband- 
detection  problem  must  take  into  account  the  dependency  of  the  calcu- 
lations on  both  frequency  and  bandwidth  of  the  signal  in  noise  to  be 
detected. 

The  single-frequency  calculation  is  satisfactory  when  the  task  is 
to  detect  a tone  or  a narrowband  of  acoustic  energy. 

A critical  term  in  the  sonar  equation  is  the  detection  threshold. 

A recent  report^  has  presented  an  excellent  and  a comprehensive  exami- 
nation of  many  recent  psychoacoustic  papers  concerned  with  the  aural 
detection  threshold.  It  has  also  presented  an  improved  f igurc-of-merit 
calculation.  The  basic  idea  of  that  calculation  will  be  used  herein  to 
develop  a signal-excess  calculation  that  will  obtain  more  accurate 
detection-range  predictions. 


Stallard,  J.N.  and  C.B.  Leslie,  "Psyclioacoust ics  and  Passive  Sonar 
Detection, " Naval  Ordnance  Laboratory  Technical  Report  74-27  (Sep  1974). 
This  report  has  an  extensive  bibliography  of  psyclioacoust ic  papers. 


The  calculation  procedure  to  be  developed  will  not  be  suitable 
for  either  manual  calculations  or  simple  desk-top  calculators.  It 
requires  a high-speed  computer.  This  is  not  really  much  of  a restriction 
In  this  day  and  age,  since  minicomputers  are  becoming  so  widespread 
In  usage  that  in  the  very  near  future  they  may  be  almost  commonplace. 

In  any  case,  information  and  knowledge  have  increased  so  rapidly  in 
recent  years — in  the  main,  through  use  of  the  computer — that  simple 
hand  calculations  are  less  and  less  adequate  for  present  day  needs. 

A review  of  the  main  results  of  Reference  1 will  be  given,  since 
the  method  for  calculating  detection  ranges  to  be  developed  herein  is 
based  on  that  work. 

DtTECTION  CAPABILITIES  OF  THE  HUMAN  EAR 

A large  number  of  psychoacoustic  investigations  have  been  made  of 
individual  factors  that  affect  the  human  ear-to-braln  detection  cap- 
abilities. However,  little  practical  application  has  been  rade  of 
these  studies  to  the  real-world  problem  of  aural  detection  through  a 
passive  sonar. 

In  Reference  1,  the  results  obtained  from  many  psychoacoustic 
papers  have  been  consolidated  into  a report  primarily  concerned  with 
quantification  of  a real-world  detection  threshold  for  the  human  ear 
and  application  thereof  to  passive  sonar.  The  work,  based  on  the 
modern  statistical  theory  of  signal  detectability,  combines  the  effects 
of  five  main  factors  to  yield  a new  detection  threshold  for  the  human 
ear  that  is  between  5 and  6 dB  higher  than  shown  in  earlier  works. 

Two  detection  thresholds  have  been  derived:  one  for  tonal  signals 

in  noise  and  another  for  broadband  signals  in  noise.  Detection  capa- 
bilities of  the  ear  are  different  for  the  two  kinds  of  signals. 

A distinctive  feature  of  the  "ear-brain"  detector  is  its  ability 
to  selectively  filter  signal  and  noise.  It  will  scarcli  out  the  region 
of  highest  signal-to-noise  ratio  of  a signal-plus-noise  spectrum,  and 
only  tlie  noise  in  that  selected  bandwidtti.  wi  1 1 affect  the  detection 
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process.  In  addition  to  the  numerous  papers  quoted  In  Reference  1,  a more 
recent  paper*"  supports  this  view  of  the  filtering  capabilities  of  the 
human  ear-brain.  The  critical  band  is  the  minimum  auditory  bandwidth; 
the  ear  forms  its  minimum  bandwidth  when  its  task  is  to  discover  or  detect 
a tone  in  noise. 

One  consequence  is  that  no  hardware  filter's  should  be  interposed  be- 
tween the  ears  of  a sonar  operator  and  output  cf  the  slgnal-plus-noise 
spectrum  of  the  sonar  when  the  task  of  the  operator  is  to  detect  a signal 
in  noise.  This  is  particularly  true  for  broadband  signals  and  noise.  A 
hardware  filter  could  seriously  degrade  the  ear-brain  detection  capabil i t ies ; 
it  could  partially  or  even  totally  exclude  the  most  detectable  portion  of 
a target  signal.  After  detection,  filtering  could  be  of  help  in  classify- 
ing and  trailing  a target. 

This  report  is  concerned  only  with  broadband  aural  detection.  The  new 
threshold  of  broadband  aural  detection  given  in  Reference  1 is 

nT=l . b+A(d' )-5  log  W (1) 

where  DT  is  the  detection  threshold  in  decibels 

A(d')  is  the  psychometric  function  in  decibels 
d'  is  tlie  detection  index 

W is  the  width  in  Hertz  of  that  portion  of  the  broadband  target 
spectrum  that  is  detected  through  the  selective  filtering 
ability  of  the  ear. 

Values  of  A(d')  are  given  in  Reference  1,  Figure  6. 

The  detection  index  depends  on  the  desired  probability  of  detection 
P CD ) and  probability  of  false  alarm  P(FA).  Tlie  detection  index  together 
with  associated  probabilities  P (D ) and  P(FA)  can  be  obtained  from  receiver- 
operator  curves  (ROC);  e.g..  Reference  1 Figure  3 or  Table  4-1. 

F.quation  (1)  indicates  that  as  W becomes  wider,  DT  is  lowered  and 
vice  versa.  A wide  band  of  broadband  signal  in  noise  is  more  detectable 
than  a narrow  band. 


2 

Moore,  B.C.,!.  and  D.H,  Raab,  "Intensitv  Discrimination  for  Noise  Rursts 
in  the  Presence  of  a Continuous,  Randstop  Rackground:  F.ffects  oi  Level, 

Width  of  the  Randstop,  and  Duration,"  Journal  of  the  Acoustical  Society 
of  America,  v.  57,  p.  400  (19751. 
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The  DT  of  Equation  (1)  is  based  on  the  results  of  a large  nuTnber  of 
psychoac rustic  experiments  and  is  thus  an  average  value  representing 
realistic  human-performance  capabilities.  There  are  many  sources  of 
variability  in  the  aural  DT  such  as  individual  sonar  operator  capability, 
experience,  fatigue,  health  condition,  motivation,  etc.  For  that  matter 
all  the  terms  of  the  sonar  equation  are  subject  to  various  degrees  of 
variability  or  fluctuation;  nevertheless,  tuere  is  a need  for  more  accu- 
rate detection-range  calculations  or,  at  least,  reasonable  approximations 
to  such.  The  method  to  be  developed  herein  is  intended  as  another  step 
toward  this  objective. 

A detection  prediction  must  be  calculated  for  a specific  bandwidth. 
The  selective  filtering  ability  of  the  car  presents  a problem,  since, 
normally,  there  is  no  way  of  knowing  in  advance  wh.<t.  signal  and  noise 
bandwidth  W will  be  most  detectable  to  the  ear.  The  ne':t  two  sections 
of  this  report  will  present  methods  for  circumventing  this  difficulty. 

For  this  purpose  Reference  1 defines 

DT'=DT+5  log  (W/af)-1.6+A(d' )-5  log  (af)  (2) 

where  ' applies  to  the  narrowbands  in  which  the  terms  of  the 
sonar  equation  are  measured  and  analyzed 

f is  ^he  geometric  center  of  W in  hertz 

a is  a proportional  constant  and  is  equal  to  0.23  for 

1 /3-octave-ana  1 ys i s bands,  0.023  for  1/30-octave  hands,  etc. 

af  is  the  bandwidth  of  the  analysis  filters;  for  constant 
bandwidth  filters,  the  term  "af"  would  be  the  constant 
bandwidth . 

The  previously  described  artifice  allows  detection  calculations  to 
be  made  without  there  being  prior  knowledge  of  tlie  band  of  frequencies 
causing  detection. 

In  Reference  1,  Equation  (2)  and  data  from  its  Figure  6 (d'  versus 
A(d'))  were  used  to  produce  the  family  of  DT'  curves  presented  herein  as 
Figure  1-1.  Tliose  curves  were  for  1/30-octave  bands;  DT'  for  1/3-octave 
bands  has  been  added  to  Figure  1-1  for  later  use. 
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FIGURE-OF-MERIT  CALCULATIONS 


The  usual  method  of  calculating  expected  detection  ranges  Is  by  using 
the  f igure-of-merlt  (FM)  form  of  the  sonar  equation. 

In  decibels,  and  following  the  notation  of  Reference  3,  the  passive 
f igure-of-merit  equation  is 


FM=SL-(NL-DI+DT)  (3) 

where  SL  is  signal-source  level 

NL  is  noise  level  at  the  receiving  sonar 
D1  is  directivity  index  of  the  receiving  sonar 
DT  is  detection  threshold. 

A better  measure  of  sonat  self-noise  is  made  at  the  sonar  output.  It 
is  called  the  equivalent  plane-wave,  sonar  self-noise  on  the  main  response 
axis  LE.  It  has  an  inherent  measure  of  the  sonar  directivity  character- 
istics at  the  time  of  measurement  and  does  not  depend  on  the  theoretical 
1)1.  It  takes  the  place  of  NL-Dt  in  the  sonar  equation.  (Note  that  it  is 
not  stated  that  LE=NL-D1.) 

In  Equation  (3),  FM  is  the  maximum  transmission  loss  that  the  signal 
can  suffer  and  still  be  detectable  at  the  receiver.  The  FM  is  used  to 
enter  a transmission  loss  versus  range  graph  by  which  FM  is  converted  to 
a detection  range. 

In  Reference  1 a method  was  discussed  for  developing  a graphical 
solution  of  the  FM  equation  that  was  intended  to  take  frequency  into 
account.  This  would  be  necessary  for  any  proper  treatment  of  the  detection 
problem,  since  all  terms  of  the  sonar  equation  are  functions  of  frequency, 
hy  using  the  method  discussed,  a graph  of  SL  was  constructed  as  a function 
of  frequency;  also,  as  a function  of  frequency,  a graph  was  constructed 
of  minimum  detectable  level  MUL. 

MDL=NL-DI+UT  (4) 


Urick,  R.J.,  "Prlr.clples  of  Underwater  Sound,"  McGraw-Hill,  Inc., 
New  York  (1975). 
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Two  MDL's  were  constructed:  one  for  tones  and  the  other  for  broad- 

band detections,  using  different  DT's  for  the  respective  cases. 

In  the  earlier  method  and  Reference  1,  this  term  was  called  an  equal 
detectability  curve;  however,  the  customary  term  MDL  will  be  used  herein. 

In  use,  the  MDL  curve  is  overlaid  on  the  SL  curve,  lining  up  the 
frequency-scaled  abscissas.  The  MDL  curve  is  then  moved  vertically 
until  it  touches  the  SL  curve.  In  this  position  the  difference  between 
the  decibel  levels  on  the  ordinates  of  the  two  graphs  is  FM — the  numerical 
solution  of  Equation  (3)  as  a function  of  frequency.  In  this  treatment 
of  the  broadband  detection  problem,  it  is  assumed  that  W is  always  the 
same  as  the  width  of  the  filters  used  to  measure  and  analyze  the  terms 
of  the  sonar  equation. 

For  broadband  detections,  it  was  assumed  in  the  previously  described 
method  that  the  bandwidth  for  the  detected  energy  was  always  one-third 
octave.  It  was  realized  in  Reference  1 that  the  detection  bandwidth  of 
the  ear  was  not  restricted  to  the  width  of  the  analysis  filters.  Thus  a 
correction  was  devised,  using  the  previously  defined  Equations  (2)  and 
(3). 

An  equation  similar  to  Equation  (3)  can  be  written  using  DT' 

FM'=SL-{NL-DI4DT' ) (5) 

Subtracting  Equation  (5)  from  Equation  (3) 

FM-FM'=nT'-DT  (h ) 

Using  Equations  (1)  and  (2) 


FM=FN'+3  log  (W/af) 


(71 
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would  determine  W of  the  detected  frequency  region,  and  a correction  to 
FM'  would  be  calculated  using  Equation  (7).  The  upper  and  lower  limits 
of  W can  be  taken  as  the  frequencies  at  which  the  separation  between  the 
SL'  and  MDL'  curves  is  3 dB  on  either  side  of  the  point  of  contact  of  the 
two  curves.  This  is  based  on  the  one-half  power  points  used  in  defining 

Ithe  bandwidth  of  filters. 

In  Reference  I,  f in  Equation  (7)  has  been  taken  as  the  point  of  con- 
tact between  the  SL'  and  MDL'  curves;  however,  it  is  more  appropriate  to 
calculate  f as  the  geometric  center  of  the  3-dB  separation  points  mentioned 
' previously.  Owing  to  the  various  SL's  and  MDL's  that  exist,  the  detection 

bandwidth  will  usually  be  unsymmetrical  about  the  point  of  contact  between 
SL'  and  MDL'  curves.  There  could  also  be  two  or  more  contact  points  be- 
tween the  curves  within  the  3-dB  limits.  It  is  possible  too,  that  the  two 
curves  could  be  in  contact  across  a large  frequency  region. 

SIGNAL-EXCESS  CALCULATIONS  FOR  BROADBAND  DETECTIONS 
With  the  previous  FM  calculations,  it  was  intended  to  take  frequency 
into  account,  since  all  of  the  terms  of  the  sonar  equation  are  freqency 
dependent.  The  method  is  adequate  for  tones  or  narrowbands;  however,  for 
broadband  detections,  the  FM  method  is  defective.  This  is  because 
transmission-loss  graphs  are  tor  practical  purposes  given  in  terms  of 
loss  at  single  frequencies.  This  means  that  the  calculated  broadband  FM 
value  must  be  used  as  a single  frequency  FM  in  order  to  enter  a graph  of 
single  frequency-transmission  loss  versus  range,  which  is  used  to  convert 
' the  FM  value  to  a detection  range. 

I 

It  would  be  impractical  to  attempt  to  provide  transmission-loss 
curves  for  all  the  possible  detection  bandwidths  that  can  occur.  These 
bandwldths,  in  practice,  are  usually  large;  so  no  single  frequency  can 
properly  represent  the  propagation  loss  or  any  of  the  other  frequency- 
dependent  terms  of  the  sonar  equation.  In  the  example  of  the  next  section, 
some  quantitative  effects  will  be  discussed.  In  addition,  W found  by 
the  FM  method  is  the  same  at  all  ranges;  whereas,  the  bandwidth  actually 
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varies  with  range,  owing  to  frequency-dependent  transmission  losses.  This 
affects  DT  because  when  W becomes  smaller,  DT  becomes  larger;  see  Equation 
(1).  A procedure  for  overcoming  these  difficulties  is  given  as  follows. 

In  terms  of  signal  excess  SE,  Equations  (3)  and  (5)  can  be  written 
as 


SE=tSL-TL)-(NL-DI+DT) 


(8) 


and 


SE'  = (SL-TL' )-(NL-DI+DT' ) (9) 

The  SL-TL  term  is  the  received  level  of  target  signal  at  the  sonar 
after  transmission  loss  (TL)  through  a range  (R). 

Each  term  of  Equations  (8)  and  (9)  is  a function  of  frequency,  now 
including  the  transmission  loss.  The  primes  have  the  same  significance 
as  before,  i.e.,  they  refer  to  the  narrow-measurement  bandwidths. 

The  procedure  now  is  to  construct  a set  of  curves  of  SL'-TL'  versus 
frequency  with  range  as  a parameter.  As  range  increases,  SL'-TL'  curves 
decrease  in  level,  the  high  frequencies  decreasing  at  a greater  rate  than 
the  low. 

Ihe  (SL'-TL* J curves  are  overlaid  with  the  MDL'  curve;  however,  this 
time,  the  ordinate  levels  of  the  two  curves  are  lined  up  as  well  as  the 
abscissa  frequencies.  The  SL'-TL'  curve  that  just  touches  the  MUL'  curve 
provides  the  range  that  would  be  obtained  if  signal  and  noise  were  limited 
to  a narrow-analysis  bandwidth  at  the  point  of  contact.  The  detection 
bandwidth  is  determined  in  the  same  way  that  the  EM  calculation  was  made 
and  a correction  is  computed  by  means  of  the  artifice  given  as  follows. 

The  signal  excess  is  zero  when  detection  can  just  occur  at  the 
probabilities  for  which  DT  has  been  determined. 

For  the  wide  detection  band.  Equation  (8)  yields  wltli  SE=0 


TL=SL-NL+DI-DT 


(10) 
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Equation  (9)  yields,  with  SE'-O 


TL'-=SL-NL4C1-DT' 


Subtracting  Equation  (11)  from  (10)  gives 

TL-TL'=DT'-DT 


(11) 


(12) 


Using  Equation  (2) 


TL=TL'+5  log  (W/af)  (13) 

The  final  range  would  then  be  found  by  sliding  MPL'  down  by  the 
amount  5 log  (W/af).  That  is,  the  wider  signal  bandwidth  associated 
with  TL  allows  a greater  transmission  loss  than  the  narrow  bandwidth  of 
TL'  and,  thus,  a longer  range. 

As  cautioned  by  Reference  1 in  the  FM  method,  the  analysis  bandwidths 
should  be  narrow  enough  so  that  tones  or  narrowbands  of  energy  can  be 
separated  from  broadband  energy,  and  W of  the  detected  signal  must  be 
greater  than  a critical  bandwidth.  Additionally,  the  analysis  filters 
should  be  narrow  enough  for  quantities  measured  through  them  to  be 
adequately  represented  at  a single  frequency. 

The  FM  method  of  Reference  1 for  calculating  detection  ranges,  while 
capable  of  being  carried  out  manually,  is  still  too  complex  and  time 
consuming  for  on-the-spot  application  such  as  in  a sonar  room. 

The  SE  method  of  calculating  detection  ranges  advocated  herein  is 
certainly  not  suited  for  manual  calculations.  A high-speed  computer  is 
essential.  This  need  will  be  appreciated  better  after  examination  of 
the  example  in  the  following  section. 

The  computer  restriction  is  not  the  severe  limitation  it  once  was. 

Use  of  minicomputers  is  growing  so  rapidly  that  they  may  shortly  be 
ubiquitous . 
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AN  EXAMPLE 


The  example  will  be  carried  out  for  1 /3-octave-ana lysis  bands. 

Smaller  bandwidths  would  be  preferable;  however,  the  only  data  that  are 
routinely  available  at  present  are  measured  through  1/3-octave  filters. 

The  receiving  sonar  has  NL  given  in  Figure  1-2.  These  are  actually 
ambient  noise  levels. 

A hypothetical  DI  is  given  in  Figure  1-3. 

The  1/3-octave  DT'  is  taken  from  Figure  1.  It  is  for  d'=l,5,  which 
corresponds  to  P(D)=0.5  and  P(FA)=0.07.  It  is  assumed  that  the  observa- 
tions will  be  made  by  a fairly  experienced  and  alert  sonar  operator,  who 
is  capable  of  working  at  the  rather  high  false-alarm  probability  of  7 
percent--a  low  threshold.  However,  when  the  operator  reports  "detection," 
it  should  be  at  a much  lower  P(FA). 

The  MPL's  constructed  from  the  previously  described  data  are  given 
in  Figure  1-4. 

A hypothetical  source  level  is  given  in  Figure  1-5. 

The  transmission-loss  model  assumed  is  for  a surface  duct.  Ducts 
are  fairly  common  in  the  oceans  of  the  world  during  the  winter  months. 
Other  models  would  be  used  as  appropriate. 

For  ranges  less  than  or  equal  to  a transition  range  (Fq) 

TL=20  logR  4-(a+')R+60  (14) 

where  TL  is  the  transmission  loss  in  decibels 
R is  the  range  in  kiloyards 

a is  the  absorption  coefficient  in  decibels  per  kiloyard 
i is  the  duct  leakage  coefficient  in  decibels  per  kiloyard 

For  ranges  greater  than  R^ 

TL=10  IorRq  +10  logR  +(a+-)R+hO  (15) 


11 


100 


Figure  l-!i  - Hypothetical 


Detectable  Levels,  Constructed  from  Figures 
1-3;  a=0.23,  d=1.5;  1/3-Octave  Bands 


Y IN  kH/ 


Thus  the  spreading  loss  Is  spherical  out  to  R_  and  cylindrical 

3 ^ 

beyond.  The  duct  transition  range  Is 


vV^ 


rH  /2(H-d) 


(16) 


where  H is  the  mixed-layer  depth 
d is  the  source  depth 

r - Cq/8 

where  r is  the  radius  of  curvature  of  sound  rays  in  the  duct 
Cq  Is  speed  of  sound  at  the  source  depth 
g is  sound  speed  gradient  with  depth. 

4 

In  the  assumed  Isothermal  duct,  Cq  Is  calculated  with  salinity  of 
35  pt;  g has  the  constant  value  0.017/sec. 

The  absorption  coefficient^  is 


1.776f ^ / 1 

32.768+f^  \l+32.  768/f^ 


0.6505f^f 
2 ~ 


0.02685f 


(17) 


where  a is  the  absorption  coefficient  in  decibels  per  kiloyard 

f is  the  frequency  in  kilohertz 

T is  the  temperature  in  degrees  Celsius 

^ . ,-(6T+118)/(T+273)  . ^ 

f^  • 21.9  X 10  kilohertz 

If  a is  the  absorption  coefficient  in  decibels  per  kiloyard,  and  a 

is  the  coefficient  in  decibels  per  kilometer,  then 

a =1. 0936a 
ra 


Leroy,  C.C.,  "Development  of  Simple  Equations  for  Accurate  and  More 
Realistic  Calculation  of  the  Speed  of  Sound  in  Sea  Water,"  Journal  of  the 
Acoustical  Society  of  America,  v.  46,  p.  216  (1969). 

^Hall,  H.R. , "Values  of  the  New  Acoustic  Absorption  Coefficient  of  Sea 
Water,"  Technical  Note  63,  Naval  Undersea  Warfare  Center  (Feb  1968). 
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The  duct  leakage  coefficient  is 

. 26.6f  (1.4)*^ 

[(1452+3. 5T)H]°*^ 


(18) 


where  i is  the  leakage  coefficient  in  decibeln  per  kiloyard 
n is  the  state  of  sea 
H is  the  duct  depth  in  feet 
f and  T are  as  previously  defined. 

The  previously  described  conversion  factor  is  also  used  to  convert 
duct  leakage  from  decibels  per  kiloyard  to  decibels  per  kilometer. 

The  T in  Equation  (18)  has  a negligible  effect;  however,  in  Equation 
(17)  at  more  than  2 kHz  it  has  a significant  effect. 

For  the  example: 

H = 400  ft  (121.9  1,1) 
d = 200  ft  (61  m) 

T = 60  F (15.6  C) 

C = 4954  ft/sec  (1510  m/s) 
r = 291.4  kyd  (266.5  km) 

Rq  = 3.6  kyd  (3.3  km) 
n = States  0,  1,  2,  4,  and  6 seas. 

The  low-frequency  cutoff  for  a 400-ft  duct  is  130  Hz,  so  it  will  not 
be  of  concern  here. 

The  amount  of  source  level  (Figure  1-5)  remaining  after  transmission 
through  various  ranges  is  given  in  Figures  1-6  through  1-10.  Each  figure 
is  for  a different  state  of  sea  in  accordance  with  Equation  (18).  Only 
the  pertinent  SL'-TE'  curves  are  given  in  the  figures,  i.e.,  the  curve 
needed  for  the  initial  range  approximation  and  the  one  needed  for  the 
corrected  range.  Naturally,  many  other  SL'-TL'  curves  were  computed  in 
the  course  of  finding  two  pertinent  for  each  case. 
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Baker,  W.F.,  "New  Formula  for  Calculating  Acoustic  Propagation  Loss 
in  a Surface  Duct  in  the  Sea,"  Journal  of  the  Acoustical  Society  of  America, 
V.  57,  p.  1198  (1975). 
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■Sea  in  1/3-Uctavc  Bands 


100 


Ranges 
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Figure  1-10  - Received  Signal  Levels  after  Transmission  through  Various  Ranges 

at  State  6 Sea  in  1 /3-Octave  Bands 
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The  procedure  used  to  calculate  a range  was  given  previously  In  the 
section  about  signal-excess  calculations.  The  procedure  for  one  example 
will  be  given  here. 

The  use  of  state  of  sea  (n)  herein  can  be  slightly  confusing.  For 
the  receiving  sonar,  n Indicates  the  Interf erlng-nolse  condition.  For 
the  received-signal  level,  n Indicates  the  amount  of  duct  leakage  the 
signal  has  suffered. 

On  Figure  1-7,  MDL'  for  a State  1 sea  has  been  "overlaid"  on  the 
SL'-TL'  curve,  which  Is  also  for  State  1 sea.  The  point  of  contact  Is 
found  to  be  4 kHz  on  the  32-kyd  (29.3  km)  SL'-TL'  curve.  The  lower  and 
upper  band  limits  (3-dB  down  rule)  are  found  to  be  1.3  and  5.8  kHz, 
respectively,  which  is  a of  4.5  kHz.  The  geometric  center  f is 
2.7  kHz,  which  Is  different  from  the  contact  point  of  4 kHz.  The  cor- 
rection from  Equation  (13)  Is  then 

5 log  4.5/0.23(2.7)=4.3  dB 

By  sliding  the  MDL'  curve  downward  4 dB  in  accordance  with  Equation 
(13),  the  new  contact  point  is  2 kHz  on  the  42-kyd  (38.4-km)  SL'-TL' 
curve.  This  satisfies  Equation  (10).  The  new  lower  and  upper  band 
limits  are  930  Hz  and  5.1  kHz,  respectively,  which  gives  a of  4.2  kHz. 

Since  the  detection  bandwidth  has  decreased  from  4.5  to  4.2  kHz, 

DT  has  increased,  and  TL  should  be  corrected  by  5 log  (4. 2/4. 5)=-0. 2 dB. 
The  detection  range  would  be  decreased  slightly;  however,  this  correction 
is  negligible.  For  MOL'  and  SL'-TL'  for  a State  0 sea,  the  secondary 
correction  was  not  negligible.  There  a first  approximation  range  of 
64  kyd  (58.5km)  was  found.  A correction  5 log  (W/af)=4  dB  yielded  a 
corrected  detection  range  at  92  kyd  (84.1  km).  However,  the  detection 
bandwidth  changed  sufficiently  for  a second  correction  of  5 log  (W^/W^ 
=-2  dB  to  yield  a final  detection  range  of  80  kyd  (73.2  km). 

Figure  1-1  la  gives  the  results  obtained  from  Figure  1-4  and  Figures 
1-6  through  1-10  for  the  five  signal  and  noise  conditions.  Detection 
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range  is  given  as  a function  of  frequency  with  n as  a parameter.  The 
bandwidth  and  frequencies  of  the  detected  signals  can  be  compared  di- 
rectly on  this  figure. 

The  bandwldths  of  the  detectable  signal  energies  in  Figure  1-1  la 
are  so  wide  that  no  single  frequency  is  representative  of  either  the 
frequency-dependent  transmission  loss  or  the  other  frequency-dependent 
terms  of  the  sonar  equation.  Even  for  the  comparatively  narrow  band 
(0.9  to  2.4  kHz)  at  80  kyd  (73,2  km)  for  a State  0 sea,  a single  fre- 
quency is  not  satisfactory.  Looking  at  Figure  1-5,  the  original  source 
level,  the  difference  in  level  between  0.9  and  2. A kHz  (the  band  limits) 
is  about  1 dB.  Looking  at  Figure  1-6,  the  signal  level  after  transmission 
through  80  kyd,  the  difference  in  level  between  these  two  frequencies  is 
now  13  dB.  Thus  12  dB  of  the  latter  difference  is  due  to  the  frequency- 
dependent  transmission  loss.  Actually,  an  "equivalent"  single-frequency 
transmission  loss  for  the  previously  described  bandwidth  and  detection 
range  could  be  selected;  however,  it  would  apply  only  to  the  specific 
SL,  MDL,  and  TL  as  well  as  80  kyd  and  the  bandwidth  from  0.9  to  2.4  kHz. 
Variations  of  the  frequency-dependent  terms  would  alter  the  range  and 
bandwidth  and,  thus,  the  selected  single  frequency.  Finally,  the 
range  and  the  corresponding  bandwidth  would  have  to  be  known  before 
making  the  single-frequency — detection-range  calculation,  which  cannot 
be  known  witliout  either  the  graphical  SE  solution  given  herein  or  an 
equivalent. 

At  the  opposite  extreme,  the  short-range,  6-kyd  detection  of  Figure 
1-1  la  shows  a 9-dB  difference  across  the  bandwidth  (2.1  to  6.9  kHz)  be- 
tween SL'  and  SL'-TL'  because  of  frequency-dependent  transmission  losses. 
The  difference  is  obtained  in  the  same  way  as  described  previously, 
using  the  bandwidth  for  n = 6 on  Figure  1-1  la,  SL'  of  Figure  1-5,  and 
SL'-TL'  of  Figure  1-10.  With  so  much  difference  across  the  detection 
bandwidth,  it  is  again  impossible  to  select  an  equivalent  single-frequency 
transmission  loss  without  prior  knowledge  of  the  bandwidth  and  frequencies 
involved. 
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In  the  first  case  (long  range,  80  kyd ) the  difference  in  level  be- 
tween the  end  points  of  the  comparatively  narrow  band  (1.5  kHz)  is 
magnified  by  the  long  range.  In  the  second  case  (short  range,  6 kyd) 
the  difference  is  owing  to  the  very  wide  bandwidth  (5  kHz). 

Figure  1-1  la  shows  that  as  the  interfering  noise  level  and  the  duct 
leakage  decrease  (state  of  sea  decreases),  the  detection  range  Increases. 

Reduced  noise  levels  allow  the  receiver  to  detect  a lower  level  signal. 

Decreased  duct  leakage  allows  signal  energy  to  be  transmitted  to  greater 
ranges;  however.  Equation  (18)  Indicates  that  there  will  always  be  some 
leakage  even  at  a State  0 sea,  and  the  leakage  Increases  directly  with 
frequency.  The  absorption  of  Equation  (17)  also  removes  high-frequency 
energy  from  the  received  signal.  Thus,  as  Figure  1-lla  shows,  long- 
range  detections  will  be  made  at  lower  frequencies  and  in  narrower  bands 
than  will  occur  at  short  ranges. 

Not  considered  herein  are  the  effects  of  varying  temperature  and 
duct  depth.  Graphs  made  from  Equation  (17)  indicate  that  as  temperature 
increases,  the  absorption  coefficient  decreases  at  frequencies  greater  ! 

than  2 kHz.  Temperature  has  a negligible  effect  on  duct  leakage.  | 

Equation  (18)  indicates  that  the  deeper  the  duct  depth,  the  less  the  | 

duct  leakage.  j 

i 

For  comparison,  detection  ranges  were  calculated  using  the  FM  3 

method  of  Reference  1.  By  overlaying  the  MDL's  of  Figure  1-A  onto  the  J 

SL'  of  Figure  1-5,  a bandwidth  of  6 kHz,  having  geometric  center  of  6 kHz,  I 

was  found  for  all  five  MDL's.  Then  5 log  (W/af)=3  dB.  The  five  FM  1 

! 

values  resulting  were  108,  98,  93,  87,  and  82  dB  for  noise  conditions 
(states  of  sea)  0,  1,  2,  3,  A,  and  6,  respectively.  To  convert  these 
FM  values  to  ranges,  the  transmission  loss  was  calculated  at  the  geo- 
metric center,  6 kHz;  see  Figure  1-12.  Entering  Figure  1-12  with  the 
five  values  of  FM,  five  ranges  were  obtained.  These  are  shown  in 
Figure  1-1  lb. 

Comparing  Figure  1-1  lb  with  1-lla,  it  is  seen  tliat  the  FM  method 
produces  smaller  detection  ranges  than  does  the  SE  mctliod.  This  bias 
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t-lgure  1-12  - Transmission  Toss  in  a Din  t at  h Kilohertz  for  Kxamplc  Given  in  Text 


toward  lower  ranges  in  the  FM  method  Is  caused  by  restricting  the 
frequency-dependent  transmission  loss  to  the  single  high  frequency  of 
6 kHz,  which  was  determined  from  the  1-yd  source  level;  see  Figure  1-5. 
However,  owing  to  the  frequency-dependent  transmission  losses,  the  de- 
tection band  moves  toward  lower  frequencies  as  range  increases.  The 
two  methods  converge  to  the  same  value  at  short  range  where  only  spread- 
ing loss  applies. 


DISCUSSION  AND  CONCLUSIONS 

On  one  hand,  the  broadband  FM  method  of  Reference  1 assumes  that  the 
detection  band  remains  the  same  at  all  ranges.  Since  this  band  is  de- 
termined by  overlaying  an  MDL'  curve  on  an  SL'  curve  (1-yd  level),  it  will 
include  high  frequencies  in  most  practical  cases.  The  single  frequency 
then  chosen  to  represent  this  band  in  a transmission  loss  versus  range 
graph  will  be  correspondingly  high.  Thus  calculated  detection  ranges 
are  biased  toward  low  values,  since  frequency-dependent  losses  are  greatest 
at  high  frequencies.  This  is  Illustrated  in  Figures  1-lla  and  1-Ilb. 

This  method  is  generally  valid  only  at  short  ranges  in  which  the  frequency 
dependent  losses  are  low.  Some  modification  could  occur  for  exceptional 
signals  and/or  receivers. 

On  the  other  hand,  a second  method  uses  values  of  SL'  and  MDL'  at  a 
single  medium  frequency  to  calculate  a figure  of  merit.  This  method  is 
difficult  to  compare  with  the  SE  method,  since  the  comparison  depends 
on  the  relative  spectral  shapes  of  SL'  and  MDL',  the  value  of  DT,  and  the 
single  frequency  chosen  for  the  calculation.  If  the  SL'  and  MDL'  curves 
are  parallel  with  negative  slopes,  if  circumstances  permit  a long  range 
detection,  if  the  detection  band  is  narrow  enough  and  can  be  reasonably 
represented  by  the  assumed  single  medium  frequency,  Chen  the  calculated 
detection  range  could  approximate  the  range  obtained  from  the  SE  method. 

At  lesser  ranges,  where  higher  frequencies  would  bo  detected,  the  medium- 
frequency  calculation  would  bias  the  calculated  ranges  toward  longer 
ranges.  The  actual  frequency-dependent  losses  would  be  greater  than 
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for  the  single-medium  frequency  method.  At  longer  ranges  than  the  pre- 
viously mentioned  long  range,  the  detection  band  would  move  toward 
lower  frequencies  than  the  fixed  medium  frequency  so  that  the  single- 
frequency  method  would  be  biased  toward  shorter  ranges  than  would  the 
SE  method.  If  SL  were  more  negatively  sloped  than  MDL,  these  effects 
would  be  intensified. 

For  the  example  given  earlier,  in  which  SL'  (Figure  1-5)  is  com- 
paratively flat,  and  MDL'  (Figure  1-4)  slopes  at  approximately  -8  to 
-9  dB/octave,  the  single-frequency  method  has  been  used  in  calculating 
detection  ranges  at  an  arbitrary  frequency  of  1.5  kHz  with  DT's  of  0 
and  -5.  The  results  are  given  in  Figure  1-1  Ic.  Two  of  the  resulting 
ranges  for  DT=-5  dB  are  equal  to  or  greater  than  the  corresponding 
ranges  calculated  from  the  SE  method.  The  other  eight  ranges  are  less 
than  for  the  SE  method.  This  single-frequency  method  can  be  made  to 
agree  with  any  individual  detection  range  from  the  SE  calculation  by 
adjustment  of  the  single  frequency  and/or  DT. 

There  seems  to  be  no  systematic  way  to  compare  the  single-frequency 
FM  method  with  the  broadband  SE  method.  The  many  circumstances  con- 
tained in  the  many  different  cases  would  have  to  be  considered  indi- 
vidually. Agreements  between  the  two  methods  of  calculating  detection 
ranges  would  generally  follow  a random  pattern. 

The  SE  method  developed  herein  does  not  have  the  previously  described 
shortcomings.  The  detection  band  is  allowed  to  vary  in  accordance  with 
the  combination  of  circumstances  that  influence  occurrence  of  a detection. 
One  possible  drawback  is  the  need  for  a high-speed  computer  for  on-the- 
spot  calculations.  However,  this  should  not  be  a problem  much  longer  in 
view  of  the  proliferation  of  high-speed  minicomputers. 

The  SE  method  also  provides  more  Information  about  an  expected 
detection.  Not  only  a detection  range  is  produced,  but  the  band  of 
frequencies  at  which  detection  is  expected  is  also  provided.  The  fre- 
quency information  can  be  quite  important  because  it  can  be  used  in 
modifying  either  the  signal  or  the  receiver  characteristics  so  as  to 
decrease  or  increase  detection  ranges. 
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All  methods  produce  the  same  detection-range  prediction  at  short 
range  when  only  spreading  loss  applies,  assuming  that  the  same  spreading- 
loss  model  is  used  in  all  methods  and  that  a correct  frequency  is  chosen 
for  the  single-frequency  FM  method.  However,  if  very  high  frequencies 
are  being  dealt  with,  the  frequency-dependent  losses  can  be  significant 
even  at  short  range.  At  long  range  the  frequency-dependent  losses  can  be 
significant  even  at  low  frequencies. 

The  method  of  computing  broadband  aural-  and  counter-detection  ranges 
developed  herein  is  too  complicated  for  expeditious  manual  calculations; 
however,  if  programed  for  a computer,  ranges  will  be  produced  of  greater 
accuracy  and,  perhaps  more  Important,  greater  consistency  than  can  be 
obtained  by  performing  the  usual  simple  manual  calculations  presently 
available  in  sonar  rooms.  Providing  for  this  method  to  be  an  element  in 
programing  the  general  purpose  digital  computers  that  are  incorporated  in 
new  sonar  systems  seems  worthwhile. 
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DTNSROC  ISSUES  THREE  TYPES  OF  REPORTS 

(1)  DTNSROC  REPORTS,  A FORMAL  SERIES  PUBLISHING  INFORMATION  OF 
PERMANENT  TECHNICAL  VALUE,  DESIGNATED  BY  A SERIAL  REPORT  NUMBER. 

(2)  DEPARTMENTAL  REPORTS,  A SEMIFORMAL  SERIES,  RECORDING  INFORMA 
TION  OF  A PRELIMINARY  OR  TEMPORARY  NATURE,  OR  OF  LIMITED  INTEREST  OR 
SIGNIFICANCE,  CARRYING  A DEPARTMENTAL  ALPHANUMERIC  IDENTIFICATION. 

(3)  TECHNICAL  MEMORANDA,  AN  INFORMAL  SERIES,  USUALLY  INTERNAL 
WORKING  PAPERS  OR  DIRECT  REPORTS  TO  SPONSORS.  NUMBERED  AS  TM  SERIES 
REPORTS;  NOT  FOR  GENERAL  DISTRIBUTION. 
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